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CN modifies the structure of the zero sound mode in the Fermi liquid in a magnetic field. 
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1 Waves in Fermi Liquids 



The standard theory of a Fermi liquid pQ allows us to predict a lot of interesting collective 
phenomena in quantum liquids, in particular, the existence of different modes of density 
fluctuations. 

The main tool of the analysis of Fermi liquids is the kinetic equation 

dn(p,x) | de(p,x) <9n(p,x) <9e(p, x) dn{p, x) = Q 
dt dp <9x dx dp 

Here and below n(p, x) is the density of particles in the phase space, and e(p, x) is the energy 

density. It is important to understand that the energy density in a specified mode is not 

equivalent to the dispersion law for the free mode, since we deal with interacting systems, 

wherein the whole system "contributes energy" into the given mode with momentum p and 

position x. The other crucial thing for understanding Fermi liquids is that the natural basis 

for the Hilbert space of excited states is the quasiparticle basis. The quasiparticle vacuum is 

the state with all energy levels below the Fermi surface are filled with particles. Thus while 

we deal with small excitations, the quasiparticles have the energy e = ep and momentum p 

confined to the Fermi sphere. 

The kinetic equation can be understood either as the semiclassical limit of the evolution 

equation for the density operator n by means of the Hamiltonian H 

%=>^ < 2 > 

with the classical Poisson bracket of {e,n} emerging from the commutator [n,H], or as the 
classical continuity equation for the current 

^ + ^ = 0, (3) 
dt ox 

where the current is normally given by 

f d 3 p <9n(p,x) f d 3 p ,de(p,x) 

j(x)= -y j2^ e{p > x) ^p— = J (2^ n(p ' x) ^^ =n(x)v ' (4) 

where we dropped the boundary terms under the integral sign. 

The basic assumption for the standard Fermi liquid theory is that the energy in the p 
mode "feels" the fluctuations of the rest of the liquid. Thus we can say that this model 
starts from a self-consistent density distribution of interacting quasiparticles, the interaction 
strength not necessarily being small. While the "equilibrium" energy density eo(p) remains 
function of momentum only, the energy fluctuation e(p, x) is a function of both momenta and 
coordinates 

e(p,x) = e (p) + £e(p,x), (5) 
and is related to the density fluctuations Sn, which, in turn, are defined as 

n(p,x) = n (p) + <5n(p,x), (6) 
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by means of a convolution with the "interaction kernel" /(p, p') 

Se(p, x) = J ^3 /(P, p')*n(P, x), (7) 

which, assuming the quasiparticle density localized on the Fermi sphere, becomes 

<5e(p,x) = J d 2 nF(8,6')5n{6',x), (8) 

where the dimensionless formfactor function F(9,9') (F = f) can ' 3e represented as a 
series in the spherical harmonics. Actually the interaction kernel / is a matrix that includes 
spin-spin interaction formfactor ~ aa'\ the integral sign in that case must be understood as 
a trace over the spin indices as well. We omit the spin indices from our analysis, since we 
limit ourselves to the class of interactions with / ~ 1; The equilibrium density and energy 
distributions no, eo are time- and coordinate-independent. If 5n describes a space-dependent 
fluctuation of particle density, the kinetic equation will become 

ddn <9eo d5n d5e driQ 
dt dp cbc <9x <9p 

Now, consider a flat-wave fluctuation of the particle density 

Sn = S{e-e F )u(B,^)e i ^ t - kr) . (10) 

Taking into account that 

<9n (p,x) xr \ 

— = _ VF j( e _ eF ) ; (11) 

we get the following integral equation from Q 

(w - kv)i/(0) = kv J ^-F(9, 9', 0, ^')v(8'), (12) 

where d 2 £l = sin9d9d<p. Restricting our analysis to (^-independent modes solely, choosing k as 
the Oz direction on the Fermi sphere, noticing that kv = kvF cos 9, introducing a convenient 
parameter s = we get 

/sin 9' ' dff 
v(9')F(9,9'). (13) 

Let us limit ourselves with the the zeroth and the first harmonics of F 

F = F + Fx (cos 9 cos 9' + sin 9 sin 9' cos(0 - $)) . (14) 



This will mean that a solution to ( 14 ) must be sought in the form 



u(8) = (Co + C lC o S 9) -, (15) 

s — cos 9 
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where Cq,Cx are some constants which are subject to conditions 



C = F / — ^— (C + Ci cos i 



sin f/f/i' „ , , , ^ cos V 

s — cos 9 ' 

(16) 

smtfdw „,„ „ „. cos 6 



„ ^ /' sin#d# 

d = F x I — - — cos 6(C + Cx cos ( 



s — cos 9 

This system of linear equations upon Cq,Cx must have a zero determinant to be solvable, 
which imposes a condition upon s thus providing us with a dispersion relation following from 
the equation 

i s lo g - l) (FoOFi + 3) + 3F lS 2 ) = Fx + 3. (17) 

The case of small Fq,Fx (both are dimensionless formfactors, describing interaction of quasi- 
particles; they can be directly inferred e.g. from the four-particle interaction Hamiltonian) 
may be called weakly interacting phase, whereas the case of large formfactors - strongly in- 



teracting. In the weakly interacting phase we expand the equation (17) around s = 1 and get 
the dispersion law 

2(F (Fi+3)+4F 1 +3) 

S = 1 + 2e fb(-Fl+3)+3fl , (18) 

whereas for the strongly-interacting case the expansion takes place around s = oo and we get 



s = 7!V Fo + ^ + — • (19) 

The oscillation mode we have described so far is known as the zero sound mode. It is important 
for us that this phenomenon is equally successfully realized both in the framework of the field 
theory (Landau Fermi liquid theory, explicated above) and in holography [2]. Of importance 
is also the fact that holography, advertised as the appropriate language for the strongly- 



coupled theory, indeed reproduces the dispersion law (19) for the large Fq,Fx, rather than 



the weakly-coupled (18) regime. In particular, even the quantum attenuation contribution 
(the absorptive part in the dispersion law not shown here) is reproduced form holography as 
well [2]. Let us emphasize that zero sound is not the single quasiparticle excitation but the 
collective excitation of the Fermi surface. 

Now we proceed to an interesting modification of zero sound due to an anomaly in the 
magnetic field. 

2 Anomaly and Chirality 

Recently Son and Yamamoto |3j have suggested that in theories with anomalies the expression 



for the current ( 24 ) be supplied with two extra terms 
f d 3 p 



Jx) 



(27T)* 



e(p , x) 9ra( P' x) - ( n dn ^ x) ) e( P , x)b - e(p , x) ( n x an(p ' x) 



dp \ dp J \ <9x 



(20) 
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Here the crucial element of the construction is the "dual" magnetic field strength fi in the 
momentum space. 

This extra current leads to the existence of several interesting effects in the presence of 
an external magnetic field. One of them is the chiral magnetic effect [3], which amounts 
to generation of an electric current along a magnetic field in presence of an axial chemical 
potential \xa 

This phenomenon was realized in the holographic approach to QCD in several different ways [5J 
El [7] and observed on lattice [8] . Holography models the chiral magnetic effect via a Chern- 
Simons term 



Scs 



J d^ vX ^A„F vX F pa , (22) 



24vr 2 

which generates the anomalous current. Recently the same term in the holographic action has 
been shown to be responsible for another very interesting effect, the chiral magnetic wave [9]. 
The essence of this effect is propagation of a chirality wave through a medium in a magnetic 
field B = (0, 0, B) 

n A ~ e*'M-fc*) > ( 2 3) 

The anomalous term in the current in the magnetic case has the immediate counterpart 
for the rotating Fermi fluid. This is the generic situation when the following substitution 
B — > 2mu> yields the correct anomaly for the rotating matter. In particular the chiral magnetic 
effect has the parallel chiral vortical effect. The presence of the proper term in the current 
in the rotating case has been remarked in [10] . In what follows we have in mind that the 
anomalous zero sound splitting can be equally considered for the magnetic or rotating cases. 

3 Anomalous Current 



The modification of the current (24) allows us to include the anomaly- driven dynamics into 
the kinetic equation. The realization of the anomaly as a contribution to the current via Berry 
phase in [3] has paved us the way to understanding how the modes of density oscillations in 
Fermi liquids are modified due to anomalies in a magnetic field. A geometrical argument 
for the appearance of the Wess-Zumino term for a Fermi surface due to a Berry phase was 
suggested recently by Zahed [UJ. Consider a Fermi liquid with right and left quasiparticles, 
described by their densities hr and ni- We have now two currents 



d 3 p 
(2^)3 



-e^ L (p,x) — Ul R>L — ) e i? , L (p,x)B- 



, <9n RiL (p,x; 
ei?,MP,x) S2 Ri L x 



<9x 



(24) 
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The oscillations of the left and right densities, that to be identical and free one from each other 
become now coupled via the anomaly term. This will lift the degeneracy in the dispersion law, 
yielding two different types of zero sound. This phenomenon will be called by us "anomalous 
zero sound". 

One might question the validity of such an approach given that we combine the non- 
relativistic physics of quasiparticles living on a Fermi surface, and an explicitly relativistic 
extra current by Son-Yamamoto. Indeed, since the extra current terms are derived from the 
Berry phase, they require the presence of a monopole at the origin of the momentum space. 
Son and Yamamoto have clearly related the extra term in the current with the anomaly. Since 
the anomalous contribution exists both for massive and massless particles we should keep the 
monopole at the origin of the momentum space for the massive particles as well as Son and 
Yamamoto did. 

Alternatively, one could appeal to the relativistic formulation of the Fermi liquid theory by 
Baym and Chin [12J, which has been shown to be much similar to the original Landau non- 
relativistic formalism. Then the argumentation by Son-Yamamoto based on Berry's phase 
would be directly applicable. 

The "dual" magnetic field fi is given by the condition that its flux through the Fermi 
surface is the topological charge 

— / dS ft = ±1, (25) 
2vr J 

whence we can choose for fi a hedgehog Ansatz 

n v = ±^, (26) 

where pp is the Fermi momentum, n is the unit normal vector on the Fermi surface, dS is the 
integration measure on the Fermi surface. 

The extra contributions in the kinetic equation ^xf 1 + ^gzr = that come from (fl ^y 1 ^ ) e(p, x)B 



are the following terms 



m 2 = ^\^ dn ^ xl )np^ 



(27) 



here the vertical arrows point out where the differential operator acts. It is clear from elemen- 
tary vector algebra identity k(k x n) = that the contribution of V e(p, x) ^fi x gn ^' x ^ j 
is zero since the wave of density as well as the wave of energy fluctuations propagate in the 
same direction of the wave- vector k. The density waves will be sought in the same form 
as prescribed by Landau theory for zero sound. We notice that such fluctuations, assuming 
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k || B, are precisely of the form that the chiral spiral wave predicted from holography 

Sn R>L (6, 0) = «5(e - e F )u R>L (0, 0)e^- kr ). 



For simplicity we stay in the axially symmetric sector v R) l = v r ^l{6)] certainly there will 
exist other zero sound modes corresponding to decomposition in higher harmonics according 
to the azimuthal angle <fi. Since we deal with a system "left fermions + right fermions + 
other degrees of freedom" , we cannot make a precise statement on how the fluctuations of the 
densities correspond to the energy fluctuations. Thus we model it as 

sin#W 

sme'dO' ( 2Q ) 



Se R = J (F s u R {e') + F A v L (9')) 
5e L = J {F AVR {6') + F s u L (e')) 



2 

The "interaction form factors" F$ , Fa contain essentially all the interesting dynamics of the 
problem. Here we chose the simplest case with F$.a = const, corresponding to the lowest 
Legendre polynomials in the spherical functions expansion. Defining them dimensionless as 
above, we allow ourselves to speak about a "weakly coupled system" where F$,a — > 0, and a 
"strongly coupled system" where F$ t A oo. Anticipating the next section of this work we 
should warn the reader that it is not at all evident that the "strongly coupled" Landau model 
of the Fermi liquid based on the nonrelativistic Fermi sphere picture and the collisionless 
kinetic equation is not a priori deemed to be equivalent to the holographic model of the same 
object. 

Under the definitions given above, the two contributions in questions become 

ikv R {9)5{e — 6f)Bvf 



(Vjfl(0))! - „ 2 



4p F 

(Vj fl (fl)) 2 = 1 2p2 F) J (FsM#) + F A y L {e')) — . 



(30) 



'F 

and analogously we can write down the contributions for the left current j^. In the equation 
for (Vj) x we have used an integration by parts trick in order to avoid nasty terms like v' R l{0)- 
Thus the equation (13) is transformed into 

( s _ cos B-b)v R {9) = (b + cos 6) (FsC r + FaC l ), 
(s - cos 9 - b)is L (9) = (b + cos 9) (F A C R + F S C L ), 
where a convenient dimensionless parameter 

has been introduced. Solving this system with regard to v r ^l and imposing the conditions 

j*M vlil ,.c R ,c l) -c l <33> 



we get thus a homogeneous system upon the coefficients Cr,l- Requiring its determinant to 
be zero we obtain the dispersion laws. They can conveniently be written down for the case 
of small and large coupling. For small coupling we see that the degeneracy has indeed been 
lifted 



l±b + e ^s(it"). (34) 



At zero magnetic field the Landau result is restored (obtainable from ( 18 ) taking F\ = 0, F$ = 
Fa = Fq). If we look at the eigenfunctions it can be shown that the axial zero sound is 
proportional to the magnetic field while the vector zero sound acquires the correction to the 
usual velocity. 

At large coupling we get 




2F S 



(35) 



Again, we comply with the Landau result (19) taking Fq = F$ 3> 1, Fs = Fa- Notice that in 
both cases the field B has been held arbitrary; the only perturbative expansion used was the 
expansion in small or large Fs a- We did not consider here the non-anomalous effects of the 
magnetic field on the structure of the modes, since in the linear order in the magnetic field 
the non-anomalous effect will be upon spin waves rather than zero sound waves. 

These dispersion laws are the main result of the field-theoretical part of this work. We 
interpret this situation in terms of two types of zero sound, axial and vector one propagating 
through a medium of left and right fermions with zero net chirality, which interact due to 
the presence of the B field. Below we analyze the two zero sounds (or the anomalous zero 
sound) from the holographic point of view. As we have remarked above the same is true for 
the rotating case. 



4 Anomalous zero sound in holography 
Anomalous splitting 

Recently a great deal of activity has taken place, advocating understanding quantum liquids 
(both of Fermi and non- Fermi type) from holography; for some of the results see |13} \TQ \15\ 
[TBI ITT] [181 f!9] ; a full review of the AdS / condensed state correspondence would certainly go 
beyond the scope of the present paper; for a pedagogic introduction see e.g. [20] . 

One can try to elucidate the nature of the mode we have observed from the holographic 
point of view. The normal zero sound (i.e. without the anomaly) in holography has been 
discovered by Karch, Son and Starinets |2j. It can be easily modified to account for the 
anomalous effects, resident in the Chern-Simons term. Consider a holographic model with 
the action 

5* = Sdbi + Sdbi + Scs ~ $cs (36) 
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where the (Dirac-Bom-Infeld-motivated) action for either left or right modes is 

Sdbi = ~ T J A^-det (g ij + 2Tta!Fh R ) (37) 



which in the metric 

„2 _ ' j 2 , / . 2 , „2jo2 



ds 2 = r^dxl 3 + (dr 2 + r 2 dQl) (38) 
in presence of the vector-potentials Aq ,R (v), and a constant field i 7 ^ = = -B becomes 



-jVVb / dr 



\ 




, i , v r4 + 52 ) (39) 

(here we have included the 2na' into the definitions of the fields) . The Chern-Simons part is 

N f 

Scs = -24^2 y drdSe^A^iV (40) 

Let us consider a system analogous to our model in the four-dimensional part of this work. 
We have already mentioned two types of holographically represented matter; now we must 
say something about the chemical potential. To comply with the previous sections we choose 
= I^R = Thus we can solve the classical equations of motion for the L and R modes 
independently by the Ansatz 

A%f( r ) = d (41) 

°' d y 1 y/fp + r 6 + r 2 B 2 V > 

where the integration constant d ~ When fluctuations 

A%' R (r) = A^ R {r) + a%' R (r, x, t) 
A 3 R ( r ) = a^' R (r,x,t) 

are taken into account, the Chern-Simons term would not have contributed to the fluctuation 
equations had three been a single type of fermions solely. Yet due to having both left and 
right modes, we can extract the VVA "anomalous vertex" structure and obtain the following 
modification of the equations of motion (25)a-c in [2] 

f3„„V7\ lf r ,,v 1 „„V\ , -iAj 



d z (/Vo) " — Mf + ?4 ) + ibaf = 0, 
dz {fgaV) + —( wt »3 + Q a o) ~ iba o' = °> 



- (43) 



d z (fga%') - ^(uaf + qa$) - ibaf = 0, 
dz {fgaf) + + qa£) + ibaf = 0, 



z 

where 

x) = J ^a L f(r, u, q )e~^, (44) 
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the remnant gauge fixing condition being 

ti V,At , V,A/ n (At .\ 
j ua + qa 3 = (J, (45) 

where we have switched to the variable z = - and introduced the functions 

r 



d 2 z 6 



and 



9 = V + / (47) 



and absorbed the normalization factor into b ~ B. Then following [2] we introduce two 
gauge-invariant field strengths 

3 (48) 
E A = + qciQ 

and construct a complex variable out of them 

E = E v + (49) 

now representing the dynamics in the chiral plane, thus we can finally write down the equation 
upon the complex variable E 



E»(A+V>(A ( f ' {z) ( 3 g 2 -^) 2 ) , m\ , E(z) {^f(zf - g 2 ) ibquzE'{z) _ 

U+ U ^ f( Z ) (,2 - ^ /(z) 2) + J + f{z) 2 -+ f{z) {q 2 _ w 2 /(z)2) " 

(50) 

This equation is then studied numerically. It is subject to the boundary equations 



EWU» = (5i) 
d z (zE(z))\ z ^oo = iujzE(z)\ z ^ 00 

The first of these conditions is the normal Dirichlet boundary condition on the boundary; the 
second one is the infalling wave condition on the horizon. These two conditions fix eigenvalues 
of this equation 

u =u(k) = a{B)k-i/3(B)k 2 . (52) 

We chose the Ansatz for u(k) containing a linear and a quadratic term only, which is confirmed 
by numeric calculations. It can be also verified by numerical analysis that the term linear in 
k is indeed purely real whereas the term quadratic in k is purely imaginary. The numerical 
dependencies of the phase velocity of the wave a and the diffusion coefficient /3 on the magnetic 
field are given in Fig. ([T]) and Fig. ^ respectively for both modes of the anomalous zero 
sound. One can readily see that the a dependence starts with a correct value, and then 

changes linearly for small values of the field. This is surely the Chern-Simons contributions, 
since the DBI could have contributed only in the order of ~ B 2 . The field lifts the degeneracy, 
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a 




Figure 2: Diffusion coefficient /3(B) of both modes of the anomalous zero sound as a function 
of the magnetic field B. 

and the splitting between the modes for small values of the field is linear for small fields; at 
large fields the DBI terms quadratic in B start to contribute. 

Thus we have produced in this section the dispersion relation for the chiral anomalous 
zero sound in holography generalizing analysis of [2]. Notice that the fluctuation effect that 
we have described exists independently from the static anomalous axial current that will be 
present in our setting as well [21] but has not been analyzed here as in the leading order the 
two effects are essentially decoupled. 

5 Discussion 

In this note we have discussed the modification of the zero sound excitation of the Fermi 
surface due to the anomaly in the external magnetic field. Since the vector and axial modes 
are coupled in the magnetic field we get two independent zero sound modes instead of the 
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single one. Using the nontrivial Berry phase in the momentum space we get the dispersion 
relations for both zero sound modes in the framework of the kinetic equations. 

Our study was focused on the peculiar longitudinal modes and a detailed analysis for the 
transverse modes can be taken into account as well. We have considered only the fluctua- 
tions of the components of the vector and axial currents ignoring the spin effects due to the 
fluctuations of the tensor currents. There should be interesting effect on the zero spin sound 
due to the mixing of scalar and tensor modes in the external magnetic field discussed in [24J. 
Another study case for us would be the involvement of the spin-spin interaction and influence 
of the magnetic field on interaction between the spin waves and sound waves. It would be also 
useful to discuss separately the effects of the possible Fermi points in the spin sound case. 

Our consideration has a lot in common with the chiral spiral wave [25J however in that case 
the modes in plasma have been investigated implying the high temperature while in our case 
we discuss the dense matter at very small or vanishing temperature. At the holographic side 
we have generalized the consideration of [23J taking into account the Chern-Simons term and 
have found the clear-cut manifestation of the anomalous contributions at the strong coupling. 
It would be interesting to investigate the possibility of the experimental observation of the 
magnetic or rotational anomalous zero sound. 
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